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pL^ ' We review and update some aspects of deuterium chemistry in the post- 

recombination Universe with particular emphasis on the formation and destruction of 
HD. We examine in detail the available theoretical and experimental data for the lead- 
ing reactions of deuterium chemistry and we highlight the areas where improvements 
in the determination of rate coefficients are necessary to reduce the remaining uncer- 
£S) \ tainties. We discuss the cooling properties of HD and the modifications to the standard 

cooling function introduced by the presence of the cosmological radiation field. Finally, 
' we consider the effects of deuterium chemistry on the dynamical collapse of primordial 

clouds in a simple "top-hat" scenario, and we speculate on the minimum mass a cloud 
must have in order to be able to cool in a Hubble time. 
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1. Introduction 



The formation of H2 and HD molecules in the post-recombination Universe plays a central role 
in the evolution of gas condensations. Even trace abundances of these molecules strongly affect 
the cooling properties of the primordial gas which would be otherwise an extremely poor radiator 
(cooling by Ly-a photons is ineffective at temperatures less than ~ 10 4 K). 

The chemistry of the early Universe has been investigated in several studies starting with the 
seminal paper by Lepp & Shull (1983). We mention in particular the work of Puy et al. (1993), Galli 
& Palla (1998; hereafter GP) and the comprehensive analysis of deuterium chemistry by Stancil, 
Lepp & Dalgarno (1998; hereafter SLD). The abundances of H2 and HD predicted at low redshift 
in these studies are of the order of 10 -6 and 10 -9 , respectively, depending on the cosmological 
model adopted (see Palla, Galli & Silk 1995 for the variation of the chemical abundances with the 
assumed baryon-to-photon ratio). A comparison of the abundance of HD obtained at z = 10 by 
GP (n[HD]/n[H] = 1.1 x 1(T 9 ), and SLD (n[HD]/n[H] = 1.6 x 1(T 9 ) for the same cosmological 
model (h = 0.5, f^o = 1, = 0.0367) shows satisfactory agreement, although this might not be 
completely significant since both calculations were based on a compilation of reaction rates obtained 
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basically from the same sources. To repeat the words of SLD, "further studies are needed to reduce 
the uncertainty in the HD abundance". One of the goals of this paper is to review the progress 
made in our understanding of HD chemistry in the three years subsequent to the publication of 
these studies. We will outline the improvements which have occurred in the meantime, and discuss 
the remaining uncertainties. 

From an observational standpoint, the abundance of atomic deuterium and molecular hydrogen 
has been measured in several cosmological clouds at redshift z ~ 2-3 in absorption along the line 
of sight to bright quasars. Deuterium in particular has attracted renewed attention because of the 
controversy about its abundance in high-redshift Ly-a clouds (see Hogan 1998 for an account and a 
resolution of the controversy) . Molecular hydrogen at high redshift was first detected by Levshakov 
& Varshalovich (1985) in a damped Ly-a system at z = 2.8. Since then, the presence of H2 has 
been confirmed in at least four additional systems (see e.g. Levshakov et al. 2000 and references 
therein). 

In Fig. 1 we compare observational data on D and H2 at high redshift with the corresponding 
abundances calculated with the standard model of GP that follows the homogenous expansion 
of the universe. The agreement between the theoretical and observed deuterium abundance is 
excellent. In the case of H2, one should keep in mind the sensitivity of this molecule to ambient 
conditions in Ly-a clouds, such as the stronger ultraviolet radiation field and the lower dust content 
with respect to local interstellar medium. These fetors can account for the considerable spread of 
observed abundances. Despite the complex phenomenology associated with damped Ly-a systems 
and the resulting uncertainty in the interpretation of the results, it is encouraging, to say the least, 
to witness the emergence of the observational foundations of the highly theoretical discipline of 
primordial chemistry. 

2. Chemical Reactions 

The formation of HD in the primordial gas follows two main routes, involving a deuteron 
exchange with H2: 

D + H 2 ^HD + H, (1) 

and 

D++H 2 ^HD + H+. (2) 

Being an isotopic modification of the the most fundamental three-electron interaction, namely, 
the H + H2 reaction, reaction (1) has received considerable interest. Thermal rate constants for this 
reaction have been measured by Ridley, Schulz, & LeRoy (1966), Westenberg & de Haas (1967), 
Mitchell & LeRoy (1973), and Michael & Fisher (1990) over a wide range of temperatures. The- 
oretical calculations employing statistical, semiclassical and quantal method have been performed 
by several groups, and show very good agreement with each other and with the experimental data. 
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The most recent studies are by Zhang & Miller (1989), Michael, Fisher, & Bowman (1990), Mielke 
et al. (1994), and Charutz, Last, <fe Baer (1997). These results are compared in the usual Arrhenius 
plot shown in Fig. 2 (the high-temperature experimental results of Michael & Fisher 1990 are not 
shown). In this paper we adopt the reaction rate computed by Mielke et al. (1994) with the DMBE 
surface (Varandas et al. 1987) which predicts much more accurate low-temperature kinetics than 
other surfaces. Their results agree with the values computed by Michael & Fisher (1990) (with the 
same surface) within ~ 20%. We also note that Zhang & Miller (1989) computed state-by-state 
cross sections for v J — > v'J' transitions and state-to-state rate constants in the temperature range 
200-1000 K. The agreement with the experimental cross section values at E ~ 1 eV (Phillips, 
Levene, &; Valentini 1989) is within a factor ~ 2. 

Reaction (2) represents the major source of HD in diffuse interstellar clouds (Dalgarno, Weisheit, 
& Black 1973). Its rate coefficient is almost constant and close to the Langevin value (2 x 
10~ 9 cm 3 s _1 , see Fig. 3). The reaction rate has been measured in the laboratory by Fehsenfeld et 
al. (1973, at T = 200 and 278 K) using a flowing afterglow technique, and by Henchman, Adams, & 
Smith (1981, at T = 205 and 295 K) using a variable-temperature selected ion flow tube technique 
(see also Smith, Adams, & Alge 1982). Gerlich (1982) performed accurate quantum-mechanical 
calculations of the rate coefficient, which are in excellent agreement with the experimental results 
of Henchmann et al. (1981) at T = 295 K, but less at T = 205 K (still within a factor ~ 2). There 
is a discrepancy with the results of Fehsenfeld et al. (1973) at both temperatures. 

The destruction of HD occurs mainly via the reverse reactions of (1) and (2), 



In general, the rate coefficients for the forwards and reverse chemical reactions are related by the 
standard thermodynamic expression (e.g. Berry, Rice & Ross 1980) 



where R is the universal gas constant, and AH , AS are the enthalpy and entropy changes. Thus, if 
the rate kf is knwon, the rate for the reverse reaction k T can be obtained directly from eq. (5). From 
the differences in the zero point vibrational energies of H2 and HD and in the ionization potentials 
of H and D, one obtains the enthalpy changes for reactions (1) and (2): AH°/R = —412 K and 
—462 K, respectively. The entropy change AS°/R can be calculated on statistical grounds (see e.g. 
Flower 2000) and is the same for both reactions: AS°/R = In 2 for para-H 2 , AS°/R = ln(2/3) for 
ortho-H2. These values are for reactant and product molecules in their ground states. Both the 
entropy and the enthalpy changes are modified when rovibrationally excited molecular states are 
involved (Flower 2000). 



HD + H^D + H 2 , 



(3) 



and 



HD + H+ -> D+ + H 2 . 



(4) 




(5) 
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Rather than using eq. (5), in this paper we prefer to provide independent fitting formulae for 
both the direct and reverse reactions that dominate the chemistry of deuterium in the primordial 
gas. For specific applications, the reader may adopt the reaction rate of the forward (or reverse, 
if better constrained) reaction, and compute the rate of the reverse (forward) reaction from the 
thermodynamic relation expressed by eq. (5). 

The rate coefficient of reaction (3) has been calculated by Shavitt (1959) using a semiempirical 
H3 energy surface. Only sparse laboratory data exist for this reaction. At T = 10 3 K, the experi- 
mental result of Boato et al. (1956) is within ~ 20% from the value calculated by Shavitt (1959). 
At lower temperatures (720-880 K), the theoretical rate is a factor ~ 2 lower than the experimental 
data by van Meersche (1951), but the extrapolation of the adopted H3 energy surface introduces 
significant uncertainty in the results at low temperatures (see Fig. 4). As a challenge to chemical 
physicists, we recall the words of Shavitt (1959): "it is unfortunate that for a reaction as basically 
important as the one considered here, the experimental data are so incomplete and inconclusive". 

Since reaction (4) is endothermic by 0.0398 eV (462 K), the removal of HD at low temperatures 
is reduced by a factor exp(— 462/T), and this can lead to significant enhancement of the HD/H2 
ratio (fractionation). The rate coefficient of this reaction has been measured in the laboratory by 
Henchmann et al. (1981) at T = 205 and 295 K, and the results are in good agreement with the 
values obtained by applying the principle of detailed balance to the reverse reaction (2). As in GP, 
we adopt the rate coefficient calculated by Gerlich (1982) over the temperature range 30-600 K 
(see Fig. 5). 

Finally, the relative abundance of D and D + is determined by the charge exchange reactions 

H+ + D -» H + D + (6) 

and its reverse 

H + D + -> H + + D. (7) 

The cross section for reaction (6) was computed by Matveenko (1974) and Hunter & Kuriyan (1977) 
for energies from 10 -3 to 7.5 eV. The two results differ by a factor ~ 2 at low energies for reasons 
unclear. Subsequent calculations by Hodges & Breig (1993), and, more recently, by Igarashi & 
Lin (1999) and Zhao, Igarashi, & Lin (2000), confirm the validity of the results of Hunter & 
Kuriyan (1977), and improve significantly the accuracy of the calculations around ~ 10 -3 eV. 
Good agreement was also found with the experimental measurements of Newman et al. (1982) and 
Esry et al. (2000). The cross section of reaction (7) was also calculated by Igarashi & Lin (1999) 
and Zhao et al. (2000). Unfortunately, no experimental data are available for this reaction. 

Watson (1976) estimated the rate of reactions (6) and (7) and Watson, Christensen & Deissler (1978) 
calculated the rate coefficient on the basis of the cross section obtained by Hunter & Kuriyan (1977). 
These rates have been widely adopted in studies of deuterium chemistry. Since reaction (6) has a 
threshold of 43 K, the rate coefficient for reaction (7) is usually obtained by multiplying that of 
reaction (6) by exp(43/T). The situation has been reanalyzed recently by Wolf Savin (2001) who 
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computed accurate rates for reactions (6) and (7) from the cross sections of Igarashi & Lin (1999) 
and Zhao et al. (2000). These results are compared to those of Watson et al. (1978) in Figs. 6 and 
7. 

In Table 1 we collect accurate fitting formulae (computed by us or by the authors quoted in 
the references) for the chemical reactions discussed in this section. These expressions update and 
replace the corresponding formulae given in Table 2 of GP. 

In addition to these reactions, additional contributions to the formation of HD in the early 
Universe come from the associative detachment reactions 

D + H -» HD + e (8) 

and 

D - + H -> HD + e, (9) 

whose rate coefficients however are not known, and can only be estimated from the corresponding 
H reactions (see SLD). Finally, minor contributions to the formation of HD come from the radiative 
association reaction 

H + D -» HD + hv, (10) 

whose rate coefficient was computed by Stancil & Dalgarno (1997), and from reactions involving 
less abundant deuterated species like HD + and H2D + (see SLD for details). 

Generally, the chemistry of HD in the primordial gas is dominated by the ion-neutral reactions 
(2) and (4) in a gas of low density (e.g. the primordial gas before the formation of the first 
structures), whereas the neutral-neutral reactions (1) and (3) become more important in conditions 
of high density and temperature (cloud collapse, shocked gas). 



3. Heating and Cooling 

In order to calculate the cooling properties of HD, one should know the population of all 
rovibrational levels. In steady-state, these are obtained by solving the balance equations 

xjJ2[RjJ'(T Tad ) + Cjj,(n,T g! , s )] =Y, x ARj'j{T rad ) + CVj(n,T gas )], (11) 
J' J' 

where J and J' indicate a generic couple of rovibrational levels. The collisional transition prob- 
abilities Cjji(n, T gas ) and Cj/j(n,T gas ) are obtained by multiplying the corresponding excitation 
coefficients, 7jj'(T gas ) and ^/j'j(T gas ), and the density of the colliding species. The terms Rjji and 
Rj'j are the radiative excitation and de-excitation rates, that can be expressed in terms of the 
Einstein coefficients Ajji and Bjj>, 

R f Ajj' + Bjj'u(ujj',T iad ), J' < J, 
JJ ' \Bjj,u(vjj,,T r! , d ), J'>J, { ) 
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where u(vjji, T rac j) is the energy density of the cosmic background radiation (CBR) per unit fre- 
quency at the temperature T ra( j, 

u{v. hV , T rad ) = MhV J 3 ' [exp(hvjj,/KT nd ) - l}' 1 . (13) 

In the primordial gas, collisional excitation of HD is dominated by collisions with H, and, 
to a less extent, with He. The coefficients for inelastic scattering of He-HD were computed by 
Green (1974) and Schaefer (1990) at temperatures T < 600 K, for < J < 3 and AJ = +l,+2. 
Collisional coefficients for rotational excitation of the system H-HD were usually derived by scal- 
ing the He-HD values with the square root of the ratio of the reduced masses of the two sys- 
tems, 7h-hd = (/^He-HD//UH-HD) 1/2 7Hc-HD, where (^h c -hd/Vh-hd) 1/2 = 1-51 (see e.g. Timmer- 
man 1996, GP). In the words of Timmermann (1996), "this assumption is however an educated 
guess, and data on H-HD are urgently needed". These data were recently provided by Flower & 
Roueff (1999) (for H and H2) and Roueff & Zeippen (1999) (for He), using full quantum- mechanical 
methods and updated energy potential surfaces to evaluate rovibrational excitation coefficients for 
collisions of HD with H, H2 and He. The results for H-HD collisions computed by Flower & 
Roueff (1999) in the temperature range 100 < T gas < 2000 K for v < 2 and J < 9 for a total of 30 
rovibrational levels are available from the CCP7 server http: //ccp7. dur.ac.uk/. 

The energy levels and the Einstein coefficients Ajj> of HD were calculated by Abgrall et 
al. (1982) who considered both dipole and quadrupole transitions and included a large number 
of rovibrational levels. Since the energy spacing of the rotational levels of HD is quite large, 
Ei/k = 128 K, E 2 /k = 383 K, E 3 /k = 764 K, etc., GP computed the HD cooling function with a 
simple four-level system (J = 0-3) adopting the collisional coefficients of Schaefer (1990). Flower 
et al. (2000) updated the calculations of GP adopting the collisional rate coefficients of Flower 
& Roueff (1999) and Roueff & Zeippen (1999). The HD cooling function computed by Flower et 
al. (2000) is also available from the CCP7 server. A useful approximation in the low-density limit 
is the expression 

AhdKH - 0)] = 2 110 (T gas )E 10 e~ E ^ kT ^ + (5/3) l21 (T gSLS )E 21 e~ E ^ kT ^ , (14) 

where E\o/k = 128 K, E 2 \jk = 255 K and the collisional rate coefficients jjji are given by Flower 
k Roueff (1999): 

7io(T gas ) = 4.4 x 10- 12 + 3.6 x 10- 13 (logT gas ) a77 cm 3 s" 1 , (15) 

and 

72i(T gas ) = 4.1 x 10" 12 + 2.0 x 10- 13 (logT gas ) - 92 cm 3 s" 1 . (16) 

The comparison between the HD cooling rate calculated by GP and Flower et al. (2000), shown 
in Fig. 8 is instructive. A simple four-level molecule is able to predict quite accurately the cooling 
rate in a wide range of temperatures (T gas < 2000 K) and densities (n[H] < 10 7 cm~ 3 ), but of course 
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fails badly in the high-temperature, high-density regime. However, for cosmological applications, 
it is important to keep in mind that the HD cooling function of Flower et al. (2000) has been 
computed assuming that the temperature of the CBR is much smaller than the gas temperature. 
This approximation is valid at low redshifts (where the two temperatures differ by more than a factor 
~ 10 for z < 20), but becomes increasingly inaccurate at higher redshifts. Since most cosmological 
scenarios of structure formation begin at z ~ 100, when T ra d 300 K, the level population of 
molecules is strongly affected by stimulated emission and absorption. In such conditions, molecules 
may become an effective heating source for the gas, because the rate of collisional de-excitation of 
the rovibrational levels is faster than their radiative decay (Khersonskii 1986; Puy et al. 1993). 

As an illustration of this effect, we plot in Fig. 9 the net heat transfer function (r — A)hd 
computed with the GP model for n(H) = 1 cm -3 as function of gas temperature at three selected 
redshifts. Note that for T gas > T rac j the cooling function is significantly decreased from the value 
computed with T ra d = 0, because of the radiative depopulation of excited states. The sudden drop 
signals the condition T gas = T ra d. Finally, for T gas < T ra d the function changes sign and becomes 
a net heating term for the gas. In cosmological simulations the heat transfer should be computed 
self-consistently at each redshift. 

4. Application: cloud collapse at z ~ 100 

In both cold dark matter and baryonic dark matter cosmological scenarios, the first objects 
predicted to enter in the nonlinear stage are the smallest ones. In cold dark matter models, it is 
expected that overdense regions with masses 10 5 -10 7 Mq first collapse in the redshift range 10 < 
z < 100 (see e.g. Cen, Ostriker, & Peebles 1993). The crucial question is whether molecular cooling 
allows the baryonic component to dissipate its kinetic energy and collapse on a timescale shorter 
than a Hubble time. This important question is fully addressed e.g. in the 1-D hydro dynamical 
calculations by Haiman, Thoul, & Loeb (1996) and in the 3-D numerical simulations presented 
by Abel, Anninos, & Norman (1997) and Bromm, Coppi, & Larson (1999). Here, we consider a 
particular issue related to the choice of the initial conditions for collapse calculations. 

Following Tegmark et al. (1997), we consider the growth of a "top-hat" overdensity region, 
an isolated spherical perturbation in a uniform density Universe (see e.g. Padmanabhan 1993). 
The radius of this region increases at a slower rate than the scale factor of the Universe (but still 
obeys the Friedman's equations) and, after reaching a maximum value (turnaround), recollapses 
to a point. The formation of a singularity is clearly an artifact of the simplified assumptions of 
the "top-hat" model. In a realistic situation, gas dynamical process (internal pressure, shocks) 
will eventually halt the collapse of the baryonic component at some finite value of the density. 
The resulting quasi-equilibrium configuration is a virialized "halo", and its subsequent evolution 
depends on the ability of the gas to cool on a timescale shorter than a dynamical timescale. 

We show in Fig. 10 the results of a sample run of the Tegmark et al. (1997) evolutionary 
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model, obtained with our chemical network and molecular cooling prescriptions. For the case 
shown, the cloud reaches virialization at z Y - u = 110, where the temperature is T v ; r = 2000 K. After 
virialization, the gas density is assumed to remain constant and uniform for the rest of the run, and 
we follow the cooling of the gas due to H2 and HD molecules. This is clearly a poor approximation 
since the cooling of the gas will of course induce an increase in the the density, and the post- 
virialization evolution of the cloud must be followed with a full hydrodynamical calculation (Galli 
et al., in preparation). Nevertheless, it is instructive to consider how rapidly the baryons are able to 
dissipate the thermal energy of the cloud via molecular cooling. As we see in Fig. 10, recombination 
reduces the ionization fraction of the cloud to negligible values, weakening the Compton cooling 
of the gas (since T ra( j < T gas ). The rapid rise of H2 (via the H - channel) around the 
sudden cooling of the gas from T v ; r down to a few hundred degrees, and the subsequent increase in 
the HD abundance further reduces the gas temperature down to few tens of degrees, establishing 
again thermal coupling of gas and cosmic radiation at redshift z ~ 10. 

Loosely speaking, a Hubble time corresponds to the redshift dropping by a factor 2 2 / 3 ~ 
1.6. For the particular case considered in this illustrative example, at a resdhift z = z v \ r /l.Q ~ 
70 the cloud temperature has dropped to T gas ~ 350, i.e. of a factor ~ 6 with respect to the 
virial temperature. Molecular cooling (mostly H2 in this range of redshift) therefore enables the 
cloud to cool significantly within a Hubble time, and eventually to collapse and form luminous 
objects. Cooling by HD appears to induce a further substantial reduction of the temperature 
at later times, but the subsequent evolution of the cloud can only be followed with a realistic 
calculation. However, it must be kept in mind that the mass scale of the fragments formed by 
the collapse of the cloud depends sensitively on the gas temperature (Mj ~ p _1 / 2 Tg/ s 2 ). The 
cooling induced by the formation of HD may in fact reduce the gas temperature below ~ 100 K 
and therefore enable the formation of primordial low-mass stars or even brown dwarfs, a possibility 
recently explored by Uheara & Inutsuka (2000) (see also Flower & Pineau des Forets 2001). These 
exciting results represent a radical cange of perspective in the field of primordial star formation 
and deserve further investigation. 



5. Conclusions 

We have examined the most effective chemical reactions leading to the formation/destruction 
of HD molecules in the early Universe, and presented a list of accurate reaction rate coefficients 
for primordial chemistry calculations that update those adopted in previous studies. We have 
analyzed the heating/cooling properties of HD molecules, stressing the relevance of a self-consistent 
calculation of the energy transfer rate between gas and radiation for cosmological applications. As 
an illustration, we have presented a simplified model for the evolution of density perturbations in 
the expanding Universe, following previous investigations by Haiman et al. (1996) and Tegmark 
et al. (1997), but including our comprehensive and updated chemical network. These results, 
together with recent findings by Uheara & Inutsuka (2000) and Flower & Pineau des Forets (2001) 
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underline the substantial contribution of HD to gas cooling during the collapse of primordial clouds. 
A preliminary conclusion from these studies is that HD is at least as important as H2 in determining 
the thermal balance of zero-metal clouds. 

This work is supported by the Italian Ministry for the University and for Scientific and Tech- 
nological Research (MURST) through a COFIN-2000 grant. It is a pleasure to thank Guillaumc 
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Table 1. Rate coefficients 



# 


reaction 


rate coefficient (cm 3 s 1 ) 


reference 


(1) 


D + H 2 - 


+ HD + H 


1.69 x io- 10 e^ 4680 / T+19880 °/ T2 ) 


Mielke et al. (1994) 


(2) 


D+ + H 2 


-> HD + H+ 


1(T 9 x [0.417 + 0.846 log T - 0.137(log T) 2 ] 


Gerlich (1982) 


(3) 


HD + H 


-^D + H 2 


5.25 x io- 1 ^^ 4430 / 1,4 " 173900 / 712 ) 


Shavitt (1959) 


(4) 


HD + H 4 


■ -> D+ + H 2 


1.1 x l - 9 e- 488 / T 


Gerlich (1982) 


(5) 


H+ + D - 


-> H + D+ 


2.00 X l -10 T 0.402 e -37.1/T _ 3 31 x 10 -17 T 1.48 


Wolf Savin (2001) 


(6) 


H + D+ -» H+ + D 


2.06 x i -io T o.396 e -33.o/r _ 2 03 x 10 -9 T -0.332 


Wolf Savin (2001) 
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Fig. 1. — Comparison of the abundances of D, D + and H2 in the primordial gas, relative to hydrogen, 
as function of redshift z calculated with the standard model of GP. Data points represent abundance 
measurements of H2 and D in damped Ly-a systems. The H2 data are taken from the list compiled 
by Levshakov et al. (2000) (detections, dots with errorbars; dashes: upper limits), whereas D data 
(triangles with errorbars) are from Buries Sz Tytler (1998a,b), O'Meara et al. (2001) and D'Odorico 
et al. (2001). 



-14- 



i 

CO 



• i — i 

o 



CD 

o 
o 

0) 



10- 



10" 



10" 



i — I — i — i — i — i — I — i — i — i — r 



B io- 



io- 

10" 

io- 



ctf 10- 



10" 
10" 




1 I I I I I I I I I I I I I I I I I I I I I I I I kJk 







2 3 4 5 

1000/T(K) 



6 



Fig. 2. — Rate coefficient for the reaction D + H2 — > HD + H according the calculations of Charutz 
et al. (1997) (short-dashed line), Zhang & Miller (1989) (long-dashed line), Mielke et al. (1994) 
(solid line). Experimental data are from Mitchell Sz LeRoy (1973) (filled squares), Westenberg & 
deHaas (1967) (open triangles), Ridley et al. (1966) (asterisks). 
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Fig. 3. — Rate coefficient for the reaction D + + H2 — > HD + H + according to the calculations of 
Gerlich (1982) (filled dots) and the measuruments by Fehsenfeld et al. (1973) (filled triangles) and 
Henchman et al. (1981) (empty triangles). The dashed line shows the Langevin value of the rate 
coefficient. The solid line shows our fit to Gerlich's results (see Table 1). 
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Fig. 4. — Rate coefficient for the reaction HD + H — > D + H2 according to the calculations of 
Shavitt (1959) for two values of the asymmetric stretching force constant A u {filled dots: A u = 
—0.358 x 10 5 dyne cm^ 1 ; empty dots: A u = —0.732 x 10 5 dyne cm -1 ). The experimental results 
of Boato et al. (1956) and van Meersche (1951) are shown by filled triangles and empty triangles, 
respectively. The solid line is our fit to the data of Shavitt (1959) for A u = —0.358 x 10 5 dyne/cm -1 
(see Table 1). 
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Fig. 5. — Rate coefficient for the reaction HD + H + — > D + + H2 according to the calculations of 
Gerlich (1982) (filled dots and the exprimental results of Henchman et al. (1981) at T = 205 and 
295 K (triangles with errorbars). The solid line shows our fit to Gerlich's results (see Table 1). 
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Fig. 6. — Rate coefficient for the reaction H + + D — > H + D + according to the calculations of 
Wolf Savin (2001) (filled dots) and Watson et al. (1978) (empty dots). The solid line shows the fit 
obtained by Wolf Savin (2001) to his results (see Table 1). 
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Fig. 7. — Rate coefficient for the reaction D + + H — > D + H + according to the calculations of 
Wolf Savin (2001) (filled dots) and Watson et al. (1978) (empty dots). The solid line shows the fit 
obtained by Wolf Savin (2001) to his results (see Table 1). 
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Fig. 8. — The cooling rate per HD molecule computed for n(H) = 1 to 10 8 cm 3 (solid lines: Flower 
et al. 2000; dashed lines: GP). Only collisions of HD with H have been considered. 
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Fig. 9. — The heat transfer function (T — A)hd for n(H) = 1 cm~ 3 versus gas temperature at 
selected redshifts. The solid line is computed ignoring the effects of the CBR (z = 0). When 
T gas > T rac j the heat exchange is a cooling term (solid lines). In the opposite case, the transfer 
becomes a heating source for the gas (dashed lines). 
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Fig. 10. — The time evolution of gas in a primordial cloud for z v ; r = 110, T v i r = 2000 K, h = 0.5, 
and fib = 0.06. The mass of the cloud is 1.35 x 10 6 M . The top panel shows the gas density n (in 
cm -3 ), the gas temperature T gas (in K) and the temperature of the cosmic background radiation 
T ra d (in K) as function of redshift z. The bottom panel shows the abundances of electrons, D, D + , 
H2 and HD, relative to H, for the same model. 



